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Abstract 

In a recent paper by Barton (J. Phys. A 40, 1011 (2007)), the 1-dimensional Klein- 
Gordon equation was solved analytically for the non-singular Coulomb-like potential 
yi(|x|) = — a/(|x| +a). In the present paper, these results are completely confirmed 
by a numerical formulation that also allows a solution for an alternative cutoff 
Coulomb potential V2(|a;|) = — |x| > a, and otherwise V2(|x|) = —a/a. 
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1 Introduction 



We consider tlie Klein-Gordon equation in one dimension, with an attractive 
vector potential V^(|a;|) (the time component of a 4- vector), and mass m. In 
units in which h = c = 1, the equation may be written [1] in the form 

{-D^ + m'^)^{x) = {E -Vf^{x), (1) 



where E denotes the energy, and, for bound states, ^ L'^{R)- Barton [2] has 
recently presented a detailed study of (1) by analytical methods in the case 
that \^ is a non-singular Coulomb-like potential with the form 



Vi{\x\) = -a/(|a;| +a). 



(2) 



where a is the fine-structure constant, and a is a positive cutoff parameter 
(called R by Barton). We shall adopt the approximation for a and the defini- 
tion of the parameter 5 given by 
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and 5 := 1/2 --/l/I^^Ri 5.32821 X 10" 




Barton's paper augments the classic work by Loudon [3] on the corresponding 
nonrelativistic problem. There is an extensive literature related to this physical 
system which is reviewed and discussed in Barton's paper. An interesting 
discussion was presented in an earlier paper by Dombey [4], who studied so- 
called hydrino states and the Klein-Gordon and Dirac Coulomb problems in 
three and two dimensions respectively. 

Since the potential V^(|a;|) is symmetric, the eigenfunctions of (1) can be taken 
to be even or odd. Most of the energies are associated with a set of odd-even 
doublets, which are essentially Balmer hke and may be expressed [2] for small 
a by the formula 
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where n— 1,2,3,..., and small- a approximations for e are given by 
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£(odd) ~ lama and efeven) a; 7 — ; — — ; — ; rr-. 

^ ' ^ ' [Q;/ma-F21og(n/2Q;ma)] 

The eigenvalues expressed by these analytical approximations of Barton are 
are not difficult to verify numerically. However, the principal goal of the present 
paper is to determine, not the Balmer-like energies, but the bottom of the 
spectrum of (1), that is to say, the energy of the 'anomalous' even state. 
Guided by the careful mathematical analysis of Ref. [2] which has exhibited 
the key spectral features generated by the potential (2), we have been able 
to devise a purely numerical search procedure which is then not limited to 
any particular design for the cutoff potential. We shall discuss this matter 
more fully in the next section. First we define the spectral parameters used 
by Barton: we adopt these for the present paper so that our results will be 
immediately comparable. Thus we define: 

ma olE 



The idea is this: the coupling a is fixed; we choose m and the cutoff radius 
a, and this determines s; the Klein-Gordon equation (1) is now completely 
specified and is solved for E for the chosen state; finally is found from 
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E and Eq.(4). For any other potential, also having a single cutoff parameter 
a (and obeying the scaling law discussed below), we may continue to use the 
same definitions. We observe that 

ds m , dp am^ 

— = — >Q and 777; = 3 > 0. 

da S dE 5{m^-E^)l 

These monotonicities mean that we can safely think of the pair {s, (3} as 
scaled images of the pair {a, E}. This observation is important since one 
of the interesting and inherent features of the problem we study is that the 
parametric curve {s, (3{s)}, s > 0, has two branches; we have to be sure a 
priori that this property is not merely a consequence of the definitions of the 
parameters. 

Another common feature exhibited by the Klein-Gordon equation, with either 
of the Coulomb-like potentials discussed in this paper, is to do with scaling. 
If we replace x hj x' = x/a, a > 0, re-write the equations, and then choose 
the scale a — 1/m, we can easily show that, for each fixed coupling a > 0, 
the energy E depends on the remaining parameters {a, m} according to the 
scaling law 

E ^ 
— = F{ma), 
m 

in which F may have more than one branch. A consequence of this law is that 
it is sufficient for the determination of the spectral function F to consider the 
problem with m fixed, say m = 1. This reduction in complexity is particularly 
important for a numerical approach. Thus, since a is in any case constant, 
what we need to do essentially is to determine how E depends on a, or, 
equivalently, how /3 depends on s. 



2 The P{s) curves 



We use fundamental numerical methods to study the lowest even eigenstate 
of (1). We make an initial guess for E in the range — ui — Ej^ <c E <C. Eu — vn 
and (for the even solution) set ■0(0) = 1 and ^'(0) = 0; discrete eigen- 
values, the wave function must be normalizeable and, in particular, satisfy 
lim2;^oo '^/'(a;) = (although need not be normalized). We then integrate nu- 
merically (by using a Runge-Kutta method with local error ~ h^) from a; = 0, 
with increasing x. We follow the size of ip^x) and count the nodes. Experience 
with Schrodinger and Sturm- Liouville eigenproblems indicates that, if E is too 
small, then |'?/'(a;)| becomes large without exhibiting the requisite number of 
nodes; meanwhile, if E is too large, there are too many nodes. Thus the en- 
ergy bounds {El, Eu} can be repeatedly and more narrowly re-chosen on the 
basis of the observed numerical behaviour of ■0(x) as x increases. This simple 
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logic works for the Klein-Gordon equation too, provided the energy sought is 
positive; when the energy sought is negative, too many nodes indicates instead 
that E is too small. This experimental observation allows us to find the lower 
branch of the curve. To be still more explicit on this point, for the upper 
branch we search in the patch (0, m), and, for the lower branch in the patch 
(— m, 0), and we use the appropriate relation between the node count and the 
size of E mentioned above. We note in passing that the interesting details of 
the graph for points near to /? = necessistated a distinct and possibly supe- 
rior approach that is described below. Our graphical results for the potential 
Vi studied by Barton are shown in Figure 1. We confirm the special values Sq 
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Fig. 1. Plot of /?(s) for the potential Vi(\x\) = -a/{\x\ + a). 
and Soo (corresponding to Barton's S2 and S4) satisfying 

f3{so) = and lim = -00. 

S ^Soo 

In terms of E and a, this limit corresponds to the supercritical limit E — > —m. 
We find for the potential Vi that sq = 0.99906868 and Soo = 6.1711. 
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Figure 1 shows the overall shape of the /3(s) curve but, as Barton found, a 
magnified patch near /3 = is needed to complete the presentation, and indeed 
to determine the minimum value = 0.99136 of s: we show this enlarged 
view in Figure 2. We note that Barton writes — s^, and also that s = Sm 
corresponds to a = a^, where ma^ = 5.28217 x 10~^. The earlier search 
policy based on the sign of E and described above was perfectly satisfactory 
for the large-scale features of the graph shown in Figure 1. However, in order 
to generate Figure 2, we were obliged to accepted the fact that a is a function 
of E (and not vice versa), even when E does not change sign. Thus, for each 
point on the graph, we fixed E in the patch (— m, m) and searched for the 
corresponding value of a that would yield a wave function i/) & L^. We found, 
for either sign of E, that the number of nodes increased with decreasing a; we 
are only interested, of course, in zero nodes, but the computer program must 
'know' how to adjust a on the basis of the values of -0(0;) for large x. At the 
end we converted to the {s, P} pair. This approach, which is independent of 
the sign of E, appears to be robust and indeed might well be used for the all 
the points on the graph. The shape of the graph of and the special values 
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Fig. 2. Enlarged view of (3{s) for the potential Vi(|x|) = — q;/(|x| + a). This graph 
shows the minimum value = 0.99136. 



{-So, -Soo, Sm} completely confirm the analytical approximations of Barton [2] 
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for this problem. 

We now consider an alternative V2 to the cutoff potential given by 



V2{\X\) = 



1^1, > a, 
\x\ < a. 



a 



(5) 



The corresponding graphs of f3{s) obtained by solving (1) with this potential 
are shown in Figs. (3) and (4). 
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Fig. 3. Plot of P{s) for the cut-off Coulomb potential V2(|a;| 
otherwise V2(|x|) = —a/a. 



—a/\x\, \x\ > a, 



The corresponding special values are sq — 1.9982289, Soo = 11.9777, and, 
from Figure 4, the minimum allowed value — 1.98216 (corresponding to 
ma^ = 1.05614 x 10"^). 
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Fig. 4. Enlarged view of j3{s) for the cut-off Coulomb potential 
V2(|x|) = — q;/|x|, \x\ > a, otherwise V2(|a;|) = —a/a. This graph shows the 



minimum value Sm = 1.98216. 
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3 Conclusion 



It should be emphasized that all the results in this paper have been obtained 
directly from Eq.(l) and are therefore completly independent of the formula- 
tion chosen by Barton for the mathematical analysis appropriate to Vi, the 
confirmation is all the more significant because of this independence. Since 
we obtain similar results for an alternative cutoff potential V2, the main fea- 
tures of the solution, discussed at length by Barton, are hkely to be inde- 
pendent of the particular form of the cutoff. This paper illustrates a con- 
structive interplay between the analytic and numerical work: the complicated 
two-branch relationship between the energy E and the cutoff radius a was 
straightforward to follow numerically once it's general form had been found 
analytically for one example. Without further mathematical work it is reason- 
ably safe to assume that the key features represented by the curves (3{s) are 
generic to one-parameter cutoff Coulomb potentials satisfying the scaling rule 
E/m = F{ma) : in particular, no discrete spectrum would be expected in the 
potential limit a — > 0. 
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